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The DAMA/Nal and DAMA/LIBRA annual modulation data, which may be in- 
terpreted as a signal for the existence of weakly interacting dark matter (WIMPs) in 
our galactic halo, are re-examined in light of new measurements of the local velocity 
relative to the galactic halo. In the vicinity of the Sun, the velocity of the Galactic 
disk has been estimated to be 250 km/s rather than 220 km/s [1]. Our analysis is 
performed both with and without the channeling effect included. The best fit regions 
to the DAM A data are shown to move to slightly lower WIMP masses. Compati- 
bility of DAMA data with nuU results from other experiments (CDMS, XENONIO, 
and CRESST I) is investigated given these new velocities. A small region of spin- 
independent (elastic) scattering for 7-8 GeV WIMP masses remains at 3(T. Spin- 
dependent scattering off of protons is viable for 5-15 GeV WIMP masses for direct 
detection experiments (but has been argued by others to be further constrained by 
Super-Kamiokande due to annihilation in the Sun). 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

Among the best motivated candidates for dark matter are Weakly Interacting Massive 
Particles (WIMPs). Direct searches for dark matter WIMPs seek to detect the scattering of 
WIMPs off of nuclei in a low-background detector. The discovery of an annual modulation 
by the DAMA/Nal experiment jj, confirmed by the new experiment DAMA/LIBRA [3| 
of the same collaboration, is the only positive signal seen in any dark matter search. The 
DAMA collaboration has found an annual modulation in its data cornpatible with the signal 
expected from dark matter particles bound to our Galactic Halo jj, [sj. Henceforth, we 
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shall use the terminology "DAMA" to refer to all of the results to date of the DAMA 
collaboration combined. Under the assumption of elastic scattering, the DAMA data alone 
is compatible with two possible WIMP mass ranges: one region around ~10-15 GeV 0, 0] 
(due to scatters primarily off of Na or channeled I scatters) and another around ~60-100 GeV 
(due to scatters primarily off of I). The details of these allowed regions will be discussed 
further below. However, other direct detection experiments, e.g. CDMS 0,9, CoGeNT [13 



COUPP [11], CRESST [I2|,ll3j, TEXONO [14|, and XENONIO [15|,ll6| have not found any 
signal from WIMPs. It has been difficult to reconcile a WIMP signal in DAMA with the 
other negative results for the case of canonical WIMPs with standard weak interactions 
motivated by Supersymmetry (SUSY). 

An alternative way to search for JA[IMP dark matter is via indirect detection of WIMP 
annihilation in the Sun or Earth 

13 



Kamiokande 



a. 



The most stringent indirect bounds are from Super- 
and will be discussed below. These bounds rely on additional physics. 



namely the strength of WIMP annihilation to neutrinos, and thus do not apply to all WIMPs. 

Previously many authors examined this question of compatibility between all the data 
sets. Prior to the DAMA/LIBRA data, Ref. [1,0] emphasized the low mass range 5-9 GeV; 
Ref. [Ill emphasized spin-dependent interactions; Ref. [H, [H, 24] studies deviations from 
the standard halo assurnptions; Ref. [1, 0] examined nonstandard velocity distributions (dark 
matter streams); Ref. [25| proposed inelastic WIMP/nucleon scattering. In the past year 
DAMA/LIBRA reported new experimental results, including presenting the modulation 
data in 36 separate energy bins and taking into account the possibility of ion channeling. 

Many authors studied the new DAMA data together with the null results of other exper- 
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2al30|,l3l|,l32| (although 



iments, and included the effects of ion channeling: 
other non-WIMP candidates were studied as well, such as mirror [33[, composite [3^, and 
WIMPless [111 dark matter. Non-neutralino WIMPs and Higgs exchange are examined in 
Ref. [s^]). In our previous work 26], we delineated the remaining regions of parameter space 
using the velocity information about our motion relative to the Galactic Halo available at 
that time. 

Recent work reexamined the local disk rotational velocities. They examined star-forming 
13, H ii and found that the fit to their data is significantly improved if the 



regions 

local velocity is taken to be higher than previously thought pL] • This shift in the velocity has 
the potential to change the interpretation of the DAMA modulation signal as well as signals 
from other experiments. It is the purpose of this paper to examine this effect. We restrict 
our studies to the standard isothermal model for the dark matter halo. We investigate both 
spin-independent and spin-dependent coupling of WIMP/nucleus elastic scattering, with 
and without channeling included. First, we find the new "best fit" regions for the DAMA 
data. Then, we compare with null results from other experiments via the likelihood method 
with goodness-of-fit. 

We begin by reviewing dark matter detection, including the role of the new velocity infor- 
mation, in Section ini In Section UTTl we review the data from DAMA and other experiments. 
We do not present much detail here, as we are using the same data and analysis techniques 
that we discussed in our previous paper [26*]. In Section HVl we apply our analysis techniques 
to obtain our results, which can primarily be found in our figures. 
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II. DARK MATTER DETECTION 



«>fmin 



WIMP direct detection experiments seek to measure the energy deposited when a WIMP 
interacts with a nucleus in the detector [40j. If a WIMP of mass m scatters elastically 
from a nucleus of mass M, it will deposit a recoil energy E = (/i^f^/M)(l — cos^), where 
fi = mM/ (m + M) is the reduced mass, v is the speed of the WIMP relative to the nucleus, 
and 6 is the scattering angle in the center of mass frame. The differential recoil rate per 
unit detector mass for a WIMP mass m, typically given in units of cpd/kg/keV (where cpd 
is counts per day), can be written as: 

where q = \/2ME is the nucleus recoil momentum, a{q) is the WIMP-nucleus cross- sect ion, 
p is the local WIMP density, and information about the WIMP velocity distribution is 
encoded into the mean inverse speed ri{E,t), 

v{E,t)=[ l^d\. (2) 

Here 

'If (3) 

represents the minimum WIMP velocity that can result in a recoil energy E and /(u, t) is 
the (time-dependent) distribution of WIMP velocities u relative to the detector. 

To determine the number of expected recoils for a given experiment and WIMP mass, we 
integrate Eqn. ([1]) over the nucleus recoil energy to find the recoil rate R per unit detector 
mass: 

rE2/Q 

R{t)= dEe{QE)-^a{q)v{E,t). (4) 

€-{QE) is the (energy dependent) efficiency of the experiment, due, e.g., to data cuts designed 
to reduce backgrounds. Q is the quenching factor relating the observed energy Ej_f,t (in 
some cases referred to as the electron-equivalent energy) with the actual recoil energy Ej-^f.'. 
Edet = QErec (explained in more detail below). The energy range between Ei and E2 is 
that of observed energies for some data bin of the detector (where experiments often bin 
observed recoils by energy). The quenching factor Q depends on the nuclear target and the 
characteristics of the detector. 

For detectors with multiple elements and/or isotopes, the total rate is given by: 

Rtot{t) = Y,fiR^it) (5) 

i 

where is the mass fraction and Ri is the rate (Eqn. (jlj)) for element /isotope i. 
The expected number of recoils observed by a detector is given by: 



Nrec = MdetTR (6) 



where M^et is the detector mass and T is the exposure time. 
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A. Cross-Section 



The a{q) cross-section term in Eqns. ([T]) & (j4]) is an effective cross-section for scatters 
witli a momentum excliange q. The momentum excliange dependence appears in form factors 
that arise from the finite size of the nucleus. The total scattering cross-section generally has 
contributions from spin- independent (SI) and spin-dependent (SD) couplings, with 



o"si + ctsd; 



(7) 



these two cross-sections are described below. 

Spin-independent (SI). For spin-independent WIMP interactions, we make the usual as- 
sumption [411 ^h^t the cross section a scales with the square of the nucleus atomic number 
A and is given by 

ao\FiE)\' (8) 



cr 



where ctq is the zero-momentum WIMP-nuclear cross-section and F{E) is a nuclear form 
factor, normalized to F{0) = 1. For purely scalar interactions. 



ao,si = -^[ZU+iA-Z)U]'. 
n 



(9) 



Here Z is the number of protons, A ~ Z is the number of neutrons, and /p and /n are the 
WIMP couplings to the proton and nucleon, respectively. In most instances, /n ~ /p; the 
WIMP-nucleus cross-section can then be given in terms of the WIMP-proton cross-section 
as a result of Eqn. ([2]): 

where the fip is the proton- WIMP reduced mass, and A is the atomic mass of the target 
nucleus. In this model, for a given WIMP mass, Up^si is the only free parameter. For the 



<70,SI — CTp.SI 



A' 



(10) 
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nuclear form factor we use the conventional Helmi form [41|, |4^ 

Spin- dependent (SD). The generic form for the spin-dependent WIMP-nucleus cross- 
section includes two couplings 4J], the WIMP-proton coupling ap and the WIMP-neutron 
coupling ttn, 



2J + 1 



[a^Sppiq) + apa^Spniq) + «n'S'„n(g)] . 



(11) 



Here, the quantities Op and a„ are actually the axial four-fermion WIMP-nucleon couplings 
in units of 1\plGp [H, IH, 47]. The nuclear structure functions Spp{q), S'„„(g), and Spn{q) 
are functions of the exchange momentum q and are specific to each nucleus. We take 



the structure functions for Aluminum from Ref. [48 
Ref. kol; and for Silicon and Germanium from Ref. 



for Sodium, Iodine, and Xenon from 
[501. These and additional structure 



functions may be found in the review of Ref. [51 



B. Velocity Distribution 

We assume the Standard Halo Model (SHM) [5|, in which the Galactic Halo is a simple 
non-rotating isothermal sphere. We use a Maxwellian distribution, characterized by an rms 



5 



velocity dispersion cr^, to describe the WIMP speeds, and we will allow for the distribution 
to be truncated at some escape velocity fesc, 



/(v) 



1 



3/2 



g-3v2/2f72 



for |v| < t>e 
otherwise. 



(12) 



Here 



erf (2;) — 2zexp(— 2;^)/7r 



1/2 



with z = Vesc/vo, is a normalization factor. The most probable speed, 



(13) 



(14) 



is used to generate unitless parameters such as z. For distributions without an escape 
velocity {v^sc 00), A^esc = 1- 

The WIMP component (halo or stream) often exhibits a bulk motion relative to us, so 
that _ 

/(u) = /(vobs + u), (15) 

where Vobs is the motion of the observer relative to the rest frame of the WIMP component 
described by Eqn. f|T2|) : this motion will be discussed below. For the truncated Maxwellian 
velocity distribution in Eqn. f[T^ . the mean inverse speed, Eqn. ([2]), becomes 



v{E,t) 



where 



1 

voy 



1 



2NcscVoy 
1 

2NescVoy 

0, 



erf(x+?/) — erf (x — y) — ^ 
erf(2;) — eif^x — y) — ■^{y + z — x)e~ 



for z < y, X < \y — z\ 
for z > y, X < \y — z\ 
for |y — 2;| < X < y + z 
for y + z < X 



(16) 



X = v„ 



(17) 



y = Vohs/ Vo , and z = Vcsc/ Vo ; 

recall fmin is given by Eqn. ([3]). Here, we use the common notational convention of rep- 
resenting 3- vectors in bold and the magnitude of a vector in the non-bold equivalent, e.g. 



|Vobs|- 



Due to the motion of the Earth around the Sun, Vobs is time dependent: Vobs 
We write this in terms of the Earth's velocity relative to the Sun as 

Vobs(i) =^Q + V(Q [£:iCOSu;(t - ti) + £2sinuj{t - ti)] , 



Vobs(i)- 



(18) 



where u = 27r/year, is the Sun's motion relative to the WIMP component's rest frame, 
= 29.8 km/s is the Earth's orbital speed, and ii and £2 are the directions of the Earth's 
velocity at times ti and ti + 0.25 years, respectively. With this form, we have neglected the 
ellipticity of the Earth's orbit, although the ellipticity is small and, if accounted for, would 
give only negligible changes in the results of this paper (see Refs. SEi for more detailed 
expressions). For clarity, we have used explicit velocity vectors rather than the position 
vectors ei and 62 used in Refs. (52I . 53] and elsewhere (the position vectors are more easily 
generalized to an elliptical orbit); the two bases are related by §1 = —62 and §2 = ^i- 
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In Galactic coordinates, where x is the direction to the Galactic Center, y the direction 
of disk rotation, and z the North Galactic Pole, 

§1 = (0.9931,0.1170,-0.01032), (19) 
£2 = (-0.0670,0.4927,-0.8676), (20) 

where we have taken Si and £2 to be the direction of the Earth's motion at the Spring 
equinox (March 21, or ti) and Summer solstice (June 21), respectively. 



C. New Measurement: Local Standard of Rest 

Unlike the Galactic disk (and the Sun), the halo has essentially no rotation; the motion 
of the Sun relative to this stationary halo is 

V0,SHM = VlsR + V0^pec , (21) 

where vlsr is the motion of the Local Standard of Rest and VQ^pec = (10, 13, 7) km/s is the 
Sun's peculiar velocity {i.e. the so-called Solar Motion). Previously, the value for the Local 
Standard of Rest was taken to be vlsr = (0, 220, 0) km/s. Now, however, new measurements 
have been made [1^, [13, [1^ and the speed that provides the best fit to data has been 
found to be higher: vlsr = (0,250,0) km/s. This change in velocity leads to the effects 
examined in this paper. For an isothermal halo, the most probable speed in Eqn. f|T^ is 
the same as the circular velocity, i.e., Vq = 250 km/s; thus the velocity dispersion is taken 
to be = 306 km/s. The best fit for the Sun's peculiar velocity may also shift slightly, 
but the changes due to this effect are much smaller and we will ignore them. Though the 
escape velocity from the Milky Way may also change (and is the subject of a future work), 
the uncertainties are so great that we continue to use the standard value v^sc = 650 km/s. 
We truncate the distribution in Eqn. (fT2|) at this value. 



D. Annual Modulation 

The count rate in WIMP detectors will experience an annual modulation as a result of 
the motion of the Earth around the Sun described above 0, [5|. Typically the count rate 
(Eqn. ([1])) has an approximate time dependence 

d/? 

— {E,t)^So{E) + Sm{E)cosuj{t-Q, (22) 

where is the time of year at which fobs(^) is at its maximum. Sq{E) is the average 
differential recoil rate over a year and Sm{E) is referred to as the modulation amplitude 
(which may, in fact, be negative). The above equation is a reasonable approximation for the 
SHM we are considering in this paper, but is not valid for all halo models, particularly at 
some recoil energies for dark matter streams; see Ref. j54] for a discussion. For the SHM, 

Sm{E) = - 



^(E, June 1) - ^(E, Dec 1) 



(23) 
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Experiments such as DAMA will often give the average amplitude over some interval, 

= ^ ^ ^ r dE SmiE). (24) 

^2 - J^l Jei 

The Earth's speed relative to the halo, fobs(^), is maximized around June 1. While this date 
varies with changes in the velocity distribution, the change is negligible (less than a day) for 
the values we are considering. 



E. Parameter Space 



Many of the parameters that factor into the expected recoil rates for a scattering detector 
are unknown, including the WIMP mass, four WIMP-nucleon couplings (SI and SD cou- 
plings to each of protons and neutrons), the local WIMP density, and the WIMP velocity 
distribution in the halo. In this paper, we shall fix the halo model to the SHM and take the 
local dark matter density po to be the estimated average density in the local neighborhood, 
0.3 GeV/cm^. In addition, we shall take /p = /n (equal SI couplings) so that there are 
only three independent scattering couplings; the SI coupling will be given in terms of the 
SI scattering cross-section off the proton, cxp^si. The parameter space we examine will then 
consist of the four parameters m, (Jp^si, ctp, and On. 



III. EXPERIMENTS 



A. DAMA 

The DAMA experiments remain the only direct detection experiments to observe a signal. 
We primarily focus on the annual modulation signal found by the DAMA group. Details of 
our analysis can be found in our previous paper (26| . 

Ion Channeling (IC): Here we review the experimental effect of ion channeling, which 
has been found to affect the results. In general only a fraction (known as the quenching 
factor) of the recoil energy deposited by a WIMP is transferred to electrons and is converted 
into useful signal in the DAMA detector {e.g. ionization or scintillation); the remainder is 
converted into phonons and heat and goes undetected. Hence measured energies must be 
corrected for this behavior to obtain the proper recoil energy, by dividing by this quenching 
factor. The DAMA detector is composed of Nal, with quenching factors Qng = 0.3 and 
Qi = 0.09. Yet, recently, the DAMA collaboration pointed out that some fraction of the 



nuclei may recoil along a crystal axis or plane [55|, in which cases the recoiling nucleus can 
travel relatively large distances along the "channels" between the crystal atoms (as opposed 
to the typical case, where the recoiling nucleus quickly collides with nearby atoms). Such 
recoiling nuclei lose very little energy to other atoms and to heat, giving nearly all their 
energy to electrons. For these cases the measured energy corresponds very nearly to the 
energy deposited by the WIMP and thus have quenching factor Q ~ 1. As a consequence, 
the detector is sensitive to lower mass WIMPs than previously thought. The new channeling 
studies revived the possiblity of DAMA's compatibility with other data sets, particularly at 
low masses. We shall examine the DAMA results both with and without including the IC 



effect. Details of our analysis can be found in our previous paper [26 
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Total Rate: While we mainly examine the DAMA experiment by analyzing the modula- 
tion amplitude S^,, DAMA can additionally constrain parameter space using the total rate 
{Sq in Eqn. (122!) ) as shown in Figure 1 of Ref. j^. The DAMA detectors do not strongly 
discriminate between WIMP scatters and background events, so an unknown and possibly 
large portion of the total observed events may be due to backgrounds (the backgrounds are 
presumed not to vary with time and do not contribute to the modulation). In the same man- 
ner as the various null experiments, we use the total number of events observed by DAMA 
to constrain the parameter space by excluding regions that would predict more events than 
observed. 

We will show regions excluded at the 90% level due to the total rate for DAMA. Con- 
straints will be shown with and without the channeling effect. Note that some choices of 
parameters that produce the correct modulation signal may actually predict a total rate 
larger than observed and, hence, those parameter choices are incompatible with the full 
DAMA data set. 



B. Null Results 

Numerous other experiments have searched for a dark matter signal, but all these exper- 
imental results, apart from DAMA, are consistent with no WIMP signal. 



Direct Detection: In our previous paper [26j we discussed in detail our use of data from 



many experiments. We do not repeat this information here. In this paper we restrict 
our analysis to the most constraining experiments: CDMS ^ [1], XENON 15, 16|, and 



CRESST I ^1^. While COUPP [11,] could be very powerful, insufficient information on the 
data runs is provided to reproduce their results and extend them to arbitrary couplings. 

For the various null experiments, we define constraints in parameter space at a certain 
exclusion level 1 — a, typically 90%, as the parameters for which the probability of seeing 
the experimental result is a. Parameters outside the corresponding contours would yield the 
observed result with a probability less than a. We say the parameters within those contours 
are compatible within the 1 — a level and parameters outside are excluded at the 1 — a level. 
The probabilities are determined using different statistical values, described in our previous 
paper, in different cases. 

Indirect Detection: An alternative way to search for WIMP dark matter is via indirect 
detection of WIMP annihilation in the Sun or Earth. When WIMPs pass through these 
objects [li,|l3, a small fraction may be captured, sink to the core, and annihilate with one 



another to produce a neutrino signal observable in Earth-based detectors such as Super- 
Kamiokande [2^ (Super-K), AMANDA [^,^, IceCube Q and ANTARES |59| \ Since 



the Sun is mostly made of hydrogen (which has spin), its capture of WIMPs depends pri- 
marily on the spin-dependent WIMP/nucleon cross section which can thus be constrained. 



The most stringent indirect bounds are from Super-K |20l . |31[ |. However, Super-K is not 
shown in the plots because their result cannot be readily reanalyzed for a modified halo. 

The Super-K constraints rely on the following assumptions: (1) the WIMP/anti- WIMP 
abundance is not highly asymmetric, which would suppress the annihilation rate; (2) the 



^ Other indirect detection methods search for WIMPs that annihilate in the Galactic Halo or near the 
Galactic Center where they produce neutrinos, positrons, or antiprotons that may be seen in detectors on 
the Earth [60,[61|. 



9 



annihilation cross-section is sufficiently high so that the capture rate of WIMPs in the Sun 
(via scattering off of nuclei in the Sun) is in equilibrium with the annihilation rate; and 
(3) the WIMP does not annihilate predominantly into the light quarks, which do not yield 
neutrinos in sufficient quantities and energies to be observed. For the Constrained Minimal 
Supersymmetric Standard Model, these assumptions are mainly satisfied for neutralinos in 
the parameter space of interest. In general, however, these assumptions need not be satisfied, 
so the Super-K constraints should be applied with caution. 



C. Analysis Techniques 

We use two statistical tests to interpret the DAMA modulation data. First, we find 
the "best ffi" regions in WIMP parameter space to fit the DAMA data by itself. More 
specifically, we use a likelihood ratio method with a global fit of four parameters to find 
the preferred parameters to produce the DAMA signal, and plot these in the plane of the 
WIMP/nucleon scattering cross section vs. WIMP mass. Second, we use a goodness-of- 
ffi test (g.o.f.) to indicate which parameters are compatible with the DAMA signal. See 
Ref. [62i] for a short review of statistics or Ref. [g^] for more extensive discussions. 



1. Best Fit to DAMA data (Likelihood Ratio Method) 

To determine the most likely parameters for producing the DAMA signal, we use the 
maximum likelihood method, based on the likelihood ratio 

L{Sm,k\fn-,(^p,Sl-,0.p-,<^-a) 

TTc F — ' ' — V 

where L is the likelihood function, Sm,k are the observed modulation amplitudes in each bin, 
and rh, (Tp^sb ^^p; and are the values of the parameters that maximize the likelihood for 
the observed Sm,k- The denominator of the above equation is the maximum likelihood value 
-^max- Confidence regions in the parameters are determined by 

2 In L(m, cXp^si, Op, On) > 2 In L^ax - 2 A In L, (26) 

where the value of A In L corresponds to the confidence level (C.L.) of the confidence region. 
Since each bin is normally distributed, this equation may instead by given in terms of the 

X^{m, 0-p,si, dp, fln) < Xmin + (27) 

where 

/ c _ cTh \2 

X [m, ap,si, ttp, On) = 2 ' (28) 

k 

(jfc is the uncertainty in 5*^,^, S'^^. = ^^^(m, CTp^si, Op, c^n) is the expected amplitude in a 
particular bin for the given parameters, and Xmm minimum value of x^- 

The Ax^ limit corresponding to different C.L.'s depends upon the number of parameters 
that have been minimized and, in the large data sample limit, is determined from a dis- 
tribution with the degrees of freedom (d.o.f.) equal to the number of minimized parameters. 
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For the four parameters considered here, Ax^ = 7.8 (90% C.L.), 16.3 {3a), 34.6 (Scr), and 
60.3 (7(t). 



The confidence region for DAMA as described here is, in fact, a 4-dimensional region in 
the (m, cTp^sij o-p, Qn) parameter space. We shall be showing 2-dimensional slices of this larger 
region for particular cases, such as the Op = an = slice corresponding to SI only scattering. 
This is not equivalent to fixing Op = = and determining a 2-dimensional confidence 
region by minimizing only over (m, o"p^si)- In the latter case, a confidence region will always 
be defined in a 2-dimensional plane while, in the former, there may be no slice in that plane 
if that plane represents a poor "fit" relative to the overall parameter space. 

Note, as will be discussed in Section HVt the confidence regions obtained via this method 
yield the most preferred parameters for producing the signal. This method, in and of itself, 
does not imply parameters outside the confidence regions are necessarily a bad fit to the 
data (or, conversely, that parameters inside these regions are a good fit) and one should be 
careful using these regions to compare with other experiments ^. 



We use the alternative statistic of goodness-of-fit (g.o.f.) to indicate when parameters are 
incompatible with DAMA data, as opposed to finding the best fit region as we did previously. 
This statistic allows for comparison of DAMA modulation data both with the total count 
rate from the DAMA data itself and with null results from other experiments (whereas 
it would be inappropriate to use the best fit regions for this purpose). To conservatively 
indicate parameters that are statistically compatible with the DAMA data set, this method 
indicates regions at which the frills within a given level using a simple goodness-of-fit 
(g.o.f.) test on the data. In contrast to the previous analysis method, there is no fit to the 
parameters here. The g.o.f. regions are defined as those parameters for which 



where Xgof the value at which the cumulative distribution function (CDF) for e.g. 
36 d.o.f. is equal to the desired level of compatibility. That is, for a desired compatibility 
level of 1 — a, there is a probability a that will exceed Xgof- Alternatively, we can say 
parameters outside of the region are excluded at the 1 — a level. 

To improve the ability of the g.o.f. test to exclude some parameters, we combine all DAMA 
bins over 10-20 keVee into a single bin (energies at which negligible signal is expected), 
resulting in a total of 17 data bins used for this test. For the 17 d.o.f. distribution, values 
of 24.8, 37.7, and 61.6 are excluded at the 90%, 3a, and 5a levels, respectively. 

While the g.o.f. regions defined here are, in fact, confidence regions — indicating the likely 
parameters to produce the DAMA signal in our theoretical framework — we do not use these 
regions in that manner. The determination of a confidence region assumes that the theoret- 
ical framework is correct, i.e. there exists some choice of parameters that is correct. That 
assumption may not be valid if, e.g., the standard halo model is not a reasonable approx- 
imation of the actual halo. Instead, we will more conservatively only use this g.o.f. test to 
exclude parameters outside of the corresponding regions. That is, parameters outside of the 



^ In statistical parlance, the determination of the confidence region for this method is decoupled from the 
goodness-of-fit. 



2. Goodness- of- Fit 




(29) 
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FIG. 1: DAMA best fit parameters (cross section and WIMP mass) for SI only scattering. The 
DAMA contom's are sliown at 5a, 3a, and 90% C.L. (outermost contour is 5a). Liglit/thin lines 
are for the old value of local standard of rest {vq = 220 km/s) while the dark/thick lines are the 
new value {vq = 250 km/s). The DAMA preferred regions are determined without (Panel A) and 
with (Panel B) the channeling effect. 

DAMA g.o.f. regions are incompatible with the DAMA signal at the given level within the 
context of this model. 



IV. RESULTS 

The likelihood ratio analysis of DAMA yields a 4-dimensional confidence region over the 
(m, (Jp^si, ^p, cin) parameter space. In Figures [H [3], andO we show the SI only (op = = 0), 
SD proton-only (cXp gi = = 0), and SD neutron-only (o"p,si = ctp = 0) slices of this 
confidence region, respectively. The thin/dashed curves are for the old velocity of the local 
standard of rest vq = 220 km/s while the solid curves are for the new local standard of rest 
vq = 250 km/s. In all cases, the first panel of any figure is without channeling, while the 
second panel is with channeling taken into account for DAMA. 

We also present results using the goodness-of-fit statistic to compare DAMA's positive 
annual modulation signal with its own bound on the total count rate plus null results from 
other experiments (CRESST I, CDMS, and XENONIO) to find what regions of parameter 
space remain compatible with all of these direct detection searches. The total rate from 
DAMA itself restricts the high mass end, and in some cases rules out the entire higher mass 
window (near 60 GeV). 

While we present our results in three particular cases — SI only, SD proton-only, and SD 
neutron-only — it is important to keep in mind that there is a larger region of parameter 
space that includes mixed couplings. 



A. Spin-Independent Couplings (SI) 



Figure [T] shows the best fit regions to the DAMA data for spin-independent couplings, 
both with and without channeling. One can see that the effect of the larger velocity is to 
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FIG. 2: Spin-independent WIMP/nucleon scattering cross section vs. WIMP mass: comparison of 
DAMA with other null experiments (CRESST I, XENONIO, and CDMS) using the goodness-of-fit 
statistic. Thin/dashed lines are for the old value of local standard of rest (vq = 220 km/s) while the 
thick/solid lines are the new value {vq = 250 km/s). The bound from the total events from DAMA 
is also shown (red solid line). Panel A does not include the effects of ion channeling while Panel B 
does take it into account. The DAMA contours are shown at 5a, 3a, and 90% C.L. (outermost 
contour is 5a). One can see that, with IC included, a sliver of parameter space remains viable at 
3a for 7-8 GeV WIMPs. 



shift the regions to slightly lower masses. Figure [2] uses the goodness-of-fit (g.o.f.) statistic 
to compare the DAMA regions with the null experiments, for the SI couplings. Most of the 
parameter space is incompatible, with the exception of a small sliver at 7-8 GeV WIMP 
masses that are still in agreement with all experiments to 3a. 
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FIG. 3: Same as Figured] but for SD proton-only scattering (a„ = 0). Note there are no 3a or 
90% C.L. contours at the higher WIMP mass region. 



B. Spin-Dependent Couplings (SD) 

Proton only [an = 0): Figure [3] shows the best fit regions to the DAMA data for spin- 
dependent couplings for proton only (a„ = 0), both with and without channeling. Again, the 
best fit regions move to slightly lower masses. In Figure H] one can see that all WIMP masses 
suitable for the DAMA data are compatible with all the other direct detection experiments. 
With channeling, the 90% region is for 7-12 GeV WIMP masses and the 3a region is for 
5-15 GeV WIMP masses. However, the COUPP results, once made available, can severely 
change these conclusions. 

The Super-K indirect detection results, which have not been plotted, can also change 



the conclusions. Ref. [Sjj have studied the constraints on neutralino WIMPs with different 
annihilation channels. If we use their results (which have not been done for the new velocity 
measurements), the conclusion for the case of SD proton-only couplings is that virtually 
all neutralinos are incompatible with the combination of SuperK plus all direct detection 
experiments; this statement assumes that there is no WIMP/antiWIMP asymmetry as well 
as the standard annihilation cross section required for thermal WIMPs. [Note that Super-K 
would rule out SD proton-only couplings, not neutron-only or SI couplings.] 

The low mass region found in the SI and the SD proton-only cases persists with a combi- 
nation of these two couplings. That is, there are portions of parameter space with both SI 
couplings and SD couphngs to the proton that are compatible with all the direct detection 
experiments. 

Neutron only {ap = 0): Figure shows the best fit DAMA regions for spin- dependent 
couplings for neutron only (a^ = 0). Figure [6] shows that there is no region of WIMP masses 
which can be compatible with all experiments in this case. While this SD neutron-only case 
is incompatible, non-zero SD couplings to the neutron are not in general incompatible: cases 
with mixed couplings (with both 7^ and a„ 7^ 0) can still be found that are consistent 
with all the experimental results. 
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10° lo' 10^ 10^ 

MwiMP (GeV) 

FIG. 4: Same as Figure [2] but for SD proton-only scattering (a„ = 0). Without Super-K included 
compatiblity is viable for 5-15 GeV WIMP masses for direct detection experiments (but has been 
argued by others to be further constrained by Super-Kamiokande due to annihilation in the Sun). 



V. CONCLUSION 

Our general result is that the new determination of the local standard of rest moves the 
best fit WIMP masses to the DAMA data to lower values. We had hoped that the DAMA 
data, which is very dependent on the WIMP velocity behavior, would change enough to allow 
new regions of parameter space compatible with all experiments. However, the bounds from 
the other experiments shifted by about the same amount, so that very little changed due to 
the analysis in this paper. 
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FIG. 5: Same as Figured] but for SD neutron-only scattering (op = 0). 
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FIG. 6: Same as Figure [2] but for SD neutron-only scattering (op = 0). No regions of parameter 
space remain. 
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